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We study the geodesic equations in the space-time of a Schwarzschild black hole pierced by an 
infinitely thin cosmic string and give the complete set of analytical solutions of these equations for 
massive and massless particles, respectively. The solutions of the geodesic equations can be classified 
according to the particle's energy and angular momentum, the ratio between the component of the 
angular momentum aligned with the axis of the string and the total angular momentum, the deficit 
angle of the space-time and as well the horizon radius (or mass) of the black hole. For bound orbits 
of massive test particles we calculate the perihelion shift, we discuss light deflection and comment 
on the Newtonian limit. 
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I. INTRODUCTION 

The motion of test particles (both massive and massless) provides the only experimentally feasible way to study the 
gravitational fields of objects such as black holes. Predictions about observable effects (light deflection, gravitational 
<sj . time-delay, the perihelion shift and the Lense-Thirring effect) can be made and compared with observations. Geodesies 
5^ ' in black hole space-times in 4 dimensional Schwarzschild space-time [l| and Kerr and Kerr-Newman space-time || 
have been discussed extensively. This has been extended to the cases of Schwarzschild-de-Sitter space-times Q as 
well as to spherically symmetric higher dimensional space-times [J] . Recently also the general solution to the geodesic 
, equation in 4 dimensional Kerr-de-Sitter Q and even general Plebanski-Demianski space-times without acceleration 
J> ■ has been found |6|. 

Cosmic strings have gained a lot of renewed interest over the past years due to their possible connection to string 
OA ■ theory 0]. These are topological defects [1] that could have formed in one of the numerous phase transitions in the 
early universe due to the Kibble mechanism. Inflationary models resulting from string theory (e.g. brane inflation) 
predict the formation of cosmic string networks at the end of inflation Q . 

Different space-times containing cosmic strings have been discussed in the past. This study has mainly been 
motivated by the pioneering work of Bach and Weyl [l(]| describing a pair of black holes held apart by an infinitely 
thin strut. This solution has later been reinterpreted in terms of cosmic strings describing a pair of black holes held 
apart by two cosmic strings extending to infinity in opposite direction. Consequently, a cosmic string piercing a 
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Schwarzschild black hole has also been discussed, both in the thin string limit [ll| - where an analytic solution can be 
given - as well as using the full U(l) Abelian-Higgs model [HGJ], 

where only numerical solutions are available. In 
the latter case, these solutions have been interpreted to represent black hole solutions with long range "hair" and are 



X 
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thus counterexamples to the No hair conjecture which states that black holes are uniquely characterized by their mass, 
charge and angular momentum. Interestingly, the solution found in [IlT | is a Schwarzschild solution which however 
differs from the standard spherically symmetric case by the replacement of the angular variable ip by (3<f>, where the 
parameter f3 is related to the deficit angle by A = 2tt(1 — /3). In this sense, the space-time is thus not uniquely 
determined by the mass, but is described by the mass and deficit angle parameter j3. 

Schwarzschild black holes pierced by cosmic strings could have formed in phase transitions in the early universe. 
One possibility would be that the Coulomb field of a charged static, spherically symmetric black hole becomes confined 
within a flux tube during such a phase transition. Interestingly, this space-time has also been used to describe the 
exterior space-time of the sun taking into account departures from perfect spherical symmetry (T3 |. 
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In order to understand details of gravitational fields of massive objects and to be able to predict observational 
consequences, it is important to understand how test particles move in these space-times. 

Geodesies in the space-time of a Schwarzschild black hole pierced by a straight cosmic string have first been 
investigated in [l5j . Spherical geodesies, i.e. geodesies in 3 spatial dimensions with constant radius r as well as 
perturbations about spherical orbits in this space-time have been discussed in [l6j and it has been observed that the 
an gula r momentum precesses around the symmetry axis of the cosmic string. The geodesies have also been discussed 
in [17| using the Hamilton- Jacobi formalism. However, neither 16] nor jl7| provide a systematic study of all possible 
geodesies. The crucial point about the solutions of the geodesic equations is that elliptic integrals have to be solved 
which was neither attempted in [l6[ nor (l7j . 

The aim of this paper is to determine the complete set of analytic solutions of the geodesic equations in the space- 
time of a Schwarzschild black hole pierced by a cosmic string and to derive analytical expressions for observable effects 
which can be used for astrophysical searches for such cosmic strings. 

Our paper is organized as follows: in Section II, we give the geodesic equations and discuss the effective potential. 
In Section III, we classify the solutions and give examples for each class. In particularly we discuss the effect of the 
conical deficit on the geodesies. Section IV contains a short discussion about the Newtonian limit, while we conclude 
in Section V. 



We consider the geodesic equation 



II. THE GEODESIC EQUATIONS 



d 2 x p rM dx p dx a _ Q 
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where T p denotes the Christoffel symbol given by 
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and s is an affine parameter such that for time-like geodesies ds 2 = g pv dx p 'dx v corresponds to proper time. The 
explicit form of the metric that we are studying in this paper is the metric of a Schwarzschild black hole pierced by a 
cosmic string [ll[ : 



ds 1 



Y,dt 2 - E^dr 2 - r 2 (d9 2 + (3 2 sin 2 



with 



1 



2 m 
r 



(3) 



This metric describes a spherically symmetric static space-time with a conical deficit angle given by 2-7r(l — /3) and 
an event horizon at the Schwarzschild radius r s — 2m. The deficit angle A is directly proportional to the energy 
per unit length fj,, which itself is equal to the tension of the string: A = 2-7r(l ~ j3) — S-kG^l. In addition we have 
to = MG, where G is Newton's constant. Note that M is a parameter that is related to the physical mass Mph ys of 
the black hole by M p h ys = 0M [16| and can hence only be interpreted as the mass of the black hole if f3 = 1. Af p h ys 
is observable through the measurement of the event horizon radius by M p h ys = (P r s)/(2G). 
The Lagrangian C for a point particle in the space-time ([3]) reads : 



C = 
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where e = for massless particles and e = 1 for massive particles, respectively. 

The constants of motion are the energy E and the component of the angular momentum of the particle that is 
aligned with the axis of the string (here the z-axis) L z [l6| : 
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dt 
ds 



constant 
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r /3sin 9— = constant . 
ds 
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In addition 



\L\ 2 = L 2 = ( ^Vr 4 + -^f- = constant . (7) 
1 1 ds sin 2 6 V ' 



Here L = (L Xl L y , L z ) is the angular momentum vector, where 



L x = -r 2 sin(/3<£) ^ - (3r 2 cos(/30) cos 9 sin (9^ (8) 
as as 



and 
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L y = r 2 cos(^0) — - /3r 2 sin(/3</>) cos 9 sin fl-r- . (9) 

The modulus of the angular momentum is always conserved, however the direction of the angular momentum is only 
conserved for j3 = 1, i.e. in the "pure" Schwarzschild case. 

Note that if we had used the Hamilton- Jacobi formalism here, the integration constant appearing in the separation 
of the equations -the so-called Carter constant- would have been equal to L 2 — L 2 . 

From variation of Q and using the constants of motion, we obtain the geodesic equations : 
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t 2 = E z [ 1 - -I , (10) 



L 2 _ L 2 

r 
L 2 



r 4 r 4 sin 2 6 ' 



(12) 
(13) 



/3 2 r 4 sin 4 6< ' 

where here and in the following the dot denotes the derivative with respect to the afhne parameter s. 

A. Effective potentials 
Equations (HU) and (12} can be rewritten in terms of effective potentials as follows: 

\r 2 + V cS (r) = , r 4 e 2 + V cS (6)=L 2 , (14) 

where 

rr , \ m L 2 L 2 m, , n . L 2 Z . . 

V cS (r) = -e- + — s , V eS (6) = . 15 

r 2r z r A sin 6 

The effective potential V e g(r) is exactly the same as that in the "pure" Schwarzschild case. Since L 2 > V e s(6), the 6 
motion is restricted: 

arcsin — ® — n — arcsin • 

Apparently, for L z = L, the motion occurs in the equatorial plane, i.e. 9 — ir/2. 

B. r(ff) and r(<fr) motion 

In the following, we will be mainly interested in the radial motion, i.e. we will solve the geodesic equations for r(9) 
and r(cj>). In order to do that we first have to eliminate the angular variable 9 from (fTB")) . By dividing (fT2"j) by ([TB^l we 

find m - 
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Psm9\/k 2 sm 2 9- 1 , (17) 
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where k 2 — {L/L z ) 2 . This can be solved to give [lq ]: 

cot 2 6» = (fc 2 - 1) sin 2 (/3<£) . (18) 



Using (HHJ), (HU), (ED and (H2J> we find: 



£ 2 - 



dS 7 / ~ L 2 sin 2 6-L 2 



r 4 sin 2 6> (19) 



and 



V<W L ? ((fc 2 - 1) sin 2 (/30) + 1 



2 



(20) 



Orbits with <j> — constant are exactly the same as in the Schwarzschild case (/3 = 1), while for orbits with (p ^ constant 
the presence of the deficit angle influences the shape of the orbits significantly. The orbits are in general non-planar 
[except for the case L z — L) and lie in a plane that has L as its normal. Since L is not conserved, this plane precesses 

^ n 

In [16| . the simplest case of non-planar orbits, namely spherical orbits and perturbations around spherical orbits 
have been discussed. Here, we want to give the complete set of possible orbits in this space-time and discuss all 
related observables. 



Classification of solutions 



Using the new variable z — =p — j and introducing 

4m 2 



we find from (HE) and (ESB using CE 
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where the third order polynomial P(z) is given by 
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This is exactly the polynomial that appears in the case of the Schwarzschild solution and the classification of solutions 
is analogue. The main difference to the "pure" Schwarzschild case is that we have non-trivial integrals on the rhs of 
(|22p. (1231) . For a given fc and B, the solutions can hence be classified according to the choice of /x and A. Obviously, 
there are only solutions for P(z) > 0. Hence, we can classify solutions according to the zeros of the characteristic 
polynomial P(z) and integrate l|22p. (f2"3")) for those z for which P(z) > 0. Note that in addition we can only integrate 
for z > —1/3 such that r > 0. Depending on the sign of the discriminant of P(z) which is given by D = g\ — 27g 2 , the 
polynomial has either three real zeros {D > 0), one real zero {D < 0) or up to two real zeros for D = which we will 
denote by a, e 2 and where e\ < e 2 < &z in the following. This is illustrated in FigfT](left) for massive test particles 
(c = 1), where the shaded regions I and II correspond to positive D, the boundary of the shaded region corresponds 
to D = and the unshaded regions III and IV correspond to negative D. Note that D > only if A < 1/3 and 
(X > 8/9. We also give the corresponding plot for massless test particles (e = 0) in Fig Q] (right). The shaded region 
corresponds to D > 0, the unshaded region to D < and the full boundary of the shaded region to D = 0. Note that 
D > (D < 0) for X/jl < A/27 (Xfi > 4/27) and that D = for A^ = 4/27 or A^ = 0. 

Comparing with (fli]). D > corresponds to the case for which UJg in < < UJg 3 *, i.e. the total "energy" E ~ £ 

of the particle is smaller than the maximum V^ ax 



and larger than the minimum V^ ln of the effective potential V c s(r). 



For D < then 
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FIG. 1: The regions corresponding to positive, negative and zero discriminant D in the A — fi plane for massive test particles 
(left) and massless test particles (right), respectively. For massive test particles the shaded regions I and II correspond to 
positive discriminant (D > 0) , the unshaded regions III and IV correspond to negative discriminant (D < 0) and the boundary 
between the two regions to zero discriminant D = 0. The vertical line at ju = 1 separates the motion of massive particles 
(e = 1) with jj, — e > and /j, — e < 0, respectively. For massless test particles, the shaded region corresponds to D > 0, the 
unshaded region to D < and the full boundary of the D > region to D = 0. 



III. SOLUTIONS OF THE GEODESIC EQUATIONS 



We can integrate ([22 



and (123 



, ( cos 6q 
arcsm 



dz 



f{4>) = --arctan[fctanO0(0- 0o))] 



dz 



such that the general solutions for the radial motion in dependence on the angular variables read 

2m 2m 

pW) - c) + i ' r(0) = p(/W>)-c) + i 



r{6) 



(25) 



(26) 



(27) 



where p is the Weierstrass elliptic function and g{9) and f(<j>) are the functions appearing in (|2"5)) and (|2"6"|) . respectively. 
Our choice of (9o depends on whether the orbit has finite maximal radius or whether it extends to infinity. For orbits 
with a finite maximal radius, we find it convenient to choose 9q = it — arcsin(l/fc) such that 



g(fi) = gi (9) = \ 



arcsm 



cos 9 



yi-fc- 2 

while for orbits that extend to infinity we choose 9q = ir/2 such that 

g(9) = g 2 {9) = i arcsin (^-^L= 
In addition we choose <fio = for all cases such that from here on 

f(4>) = — — arctan [kt&n(/3<j>)] . 



(28) 



(29) 



(30) 



The constant c in (|2~7T) is given by c = J 



— dz For D < we can always choose c = 0. For D > we choose 

c = for zq = oo, c = U3\ for zq = ei, c = ui + for zq = e 2 and c — u>2 for zq = Z3. Here 



%/ei - e 3 



CJ 2 



. K(/CQ 
' y/ei - e 3 



(31) 
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where K denotes the complete elliptic integral of the first kind and JC is the modulus of the elliptic integral with 

K = FEEH . (32) 

V ei - e 3 

Moreover we have JC' = yjl — JC 2 . For D > 0, D = and D < there are various replacements of the Weierstrass 
function p by real valued Jacobi elliptic functions [la ]. 

In the following, we will discuss different type of orbits, which we denote as follows : 

1. Bound orbit : an orbit for which r varies between two finite values. 

2. Spherical orbit : a special bound orbit for which r is constant. 

3. Bound terminating orbit : a bound orbit that ends at the singularity r = 0. 

4. Unbound terminating orbit : an orbit for which r varies between r = and infinity. 

5. Escape orbit : the orbit comes from infinity, approaches a finite r and goes again to r = oo. 

A. Geodesies for D > and E 2 - e > 

This case corresponds to the choice V^ n < — ~ e < V^ ax and E 2 — e > (region I in the A-/i-plot). The effective 

potential and the corresponding characteristic polynomial P(z) are shown in Fig|5J The line E ~ e intersects the 
effective potential V c g(r) twice. The polynomial P(z) has in fact three real zeros e\, e2, e%. However, for < z < e% 
we have P(z) < and in addition e± < —1/3. We are hence allowed to integrate from e-2, to z = —1/3, which 
corresponds to integration from r = r 2 to r = oo and from to z — oo which corresponds to integration from r = 
to r = f3. 

1. Escape orbit and light deflection 
The range of allowed values for r is indicated by the red horizontal line in FigJ5]and is r 2 < r < oo. 

Effective potential Characteristic polynomial 




FIG. 2: The effective potential V e a{r) (left) and the characteristic polynomial P(z) (right) of a massive test particle s = 1 for 
D > 0. The red horizontal line corresponds to the range of r values for an escape orbit. 

The plots of escape orbits of a massive test particle and a massless test particle are given in FigJ3] and FigJH 
respectively. Obviously, for f3 close to one, test particles get simply deflected. However, if (3 is significantly smaller 
than one, particles will approach the minimal radius r2, then whirl the black hole and finally move off to infinity 
again. This is particularly important for the deflection of light. 
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FIG. 3: We show an escape orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin#, z = rcosO) (left), in 
the plane perpendicular to L (middle) and in K 3 (right), (a), (b) and (c) correspond to /3 = 0.33, /3 = 0.5 and (3 — 0.99, 
respectively. Here, we have chosen E = 1.2, L z — 5, k = 1.2 and m = 1. The yellow circles in the /9-z-plane and in the plane 
perpendicular to L correspond to circles with radius n, respectively. 



In order to see the effect of (3 on the deflection of light (e = 0) we set k — 1 (i.e. we choose 9 — ir/2) and select the 
initial point at rp = r%. The requirement r(cj>) = oo imposes a bound on ip which reads: 



1 



V e l - e 3 Jo 
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l-K, 2 sin 2 (<^) 



(33) 
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(a) 



with = 0.33 





(b) 



with P = 0.50 





30 20 




(c) 



with P = 0.96 



FIG. 4: We show an escape orbit for a massless test particle (e = 0) in the p-z-plane (with p = rsinO, z — rcos9) (left), 
in the plane perpendicular to L (middle) and in R 3 (right), (a), (b) and (c) correspond to ft = 0.33, /3 = 0.5 and ft — 0.96, 
respectively. Here, we have chosen E = 0.75, L z = 5, k = 1.1 and m = 1. The yellow circles in the /9-z-plane and in the plane 
perpendicular to L correspond to circles with radius ri, respectively. 



Therefore the deflection angle A<fi is given by: 



dip 



2wi 



■xl,-l 



Vei~e 3 Jo y J 1 _ !C 2 sin 2 (¥ ,) 
The term in the square bracket of is the expression of A(f> for the "pure" Schwarzschild case, i.e. 



(34) 
1 and the 
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second term takes care of the fact that the space-time has a conical deficit. 

Using the parametrized Post- Newtonian (PPN) formalism |20(, which describes deviations from standard General 
Relativity, the angle of light deflection can be given by A<f> = |(1 + 71) 1. 75" = \{1 + ji)(A<j))s assuming that 
the parameter 71 is equal to unity for General Relativity. (A0)g denotes the General Relativity value assuming, the 
exterior space-time of the massive body to be given by the Schwarzschild solution. Different experimental tests [2l| - l23j 
have given a value of 71 — 1 = (2.1 ± 2.3) ■ 10~ 5 , hence A ^^^ s < 10~ 5 . If we now assume that the deviations of 
Acf> from the Schwarzschild value are not due to a modification of General Relativity, but due to the presence of a 
cosmic string, we can approximate the energy per unit length of the cosmic string. Reinstalling factors of c 2 we then 
find (1 — /3) < 10 -11 which for the deficit angle gives A < 10~ 10 . This transfers to a bound on the energy per unit 
length given by : fj, < 10 16 ^. 



2. Bound terminating orbit 

The range of allowed r values is indicated in Fig 12] and is < r < n. We don't discuss this case in detail here, 
because it is qualitatively similar to the case with D > and E 2 — e < (see discussion below) . 



B. Geodesies for D > and ^f^- < 

This case corresponds to the choice V^ ln < E ' 2 ~ e < V™ x and E 2 — e < (region II in the A- //-plot). The effective 

potential and the corresponding characteristic polynomial P{z) are shown in Fig. [5] The line — 2 ~ £ intersects the 
effective potential V Q g(r) three times at r±, and r%. These intersection points correspond to the zeros of the 
characteristic polynomial, which we denote by e\, ei and 63. Apparently, we are allowed to integrate between 00 and 
d, which corresponds to integration from r = to r = n, and from 62 to e3, which corresponds to integration from 
r = T2 to r = r%. 



Effective potential Characteristic polynomial 
0.15i . 1 . 1 . 1 1 0.05 . . . . r- 




FIG. 5: The effective potential V D s(r) (left) and the characteristic polynomial P(z) (right) of a massive test particle e = 1 for 
D > 0. The red horizontal line corresponds to the range of r values for a bound orbit, while the yellow line corresponds to the 
range of r for a bound terminating orbit. 



1. Bound terminating orbit 
The range of allowed r values is indicated by the yellow line in FigJS] and is < r < r\ . 

A plot of a bound terminating orbit for a massive test particle is given in Fig. [6] We emphasize on the deformation 
of the orbits when changing the deficit parameter (3. While for j3 = 0.99, the orbit is still nearly planar, this changes 
with the decrease of j3, i.e. the increase of the deficit angle. While the orbit looks heart-shaped for j3 = 0.99, the 
increase of the deficit angle leads to a pretzel-like structure. 
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(c) with p = 0.99 



FIG. 6: We show a bound terminating orbit for a massive test particle (s = 1) in the p-z-plane (with p — rsin#, z — rcos9) 
(left), in the x-y-plane (middle) and in R (right), (a), (b) and (c) correspond to /3 = 0.25, /3 = 0.75 and /3 = 0.99, respectively. 
Here, we have chosen E — 0.9975, L z — 5, k — 2.5 and m = 1. The blue and the dark yellow circles in the p-z-plane and in the 
£-j/-plane correspond to circles with Schwarzschild radius r s = 2m = 2and maximal radius n, respectively. 



2. Bound orbit and perihelion shift 
The range of allowed r- values is indicated by the red line in Fig jS] and is r-i < r < r^. 

A bound orbit of a massive test particle is given in Fig. [7] for different choices of (3. For /3 close to one, the orbit 
corresponds to a standard closed orbit, while decreasing (3 leads to the appearance of additional loops in the orbit 
(see the cases (3 — 0.25 and j3 — 0.75). This is important in order to understand the perihelion shift of test particles 
orbiting a massive body that is pierced by a cosmic string. 
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FIG. 7: We show a bound orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin#, z — rcosO) (left), in 
the plane perpendicular to L (middle) and in R 3 (right), (a), (b) and (c) correspond to f3 = 0.25, /3 = 0.75 and f3 = 0.99, 
respectively. Here, we have chosen E = 0.9975, L z — 5, k = 2.5 and m — 1. The yellow and magenta circles in the p-z-plane 
and in the plane perpendicular to L correspond to circles with radius ri and maximal radius rs, respectively. 



For bound orbits the perihelion shift depends on the deficit parameter f3. This effect has been calculated approx- 
imately in [lij . Here, we will give the exact analytic formula for the perihelion shift. We choose k = 1 in order to 
study planar orbits with 9 = ir/2. We then have for the motion in the plane: 

rU) = : . (35) 
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When the test particle moves from r% — > r% — > ?'3 we find for the change in 0: 



p sje\ - e 3 



such that the perihelion shift A(j) reads 



A<f> = 6<f> - 2n = 4^E=L-2vr . (37) 

P V e l - e 3 

In FigEl we plot the value of A<fi in dependence on /3 for m = 1, £ z = 5, k = 1 and _E = 0.999. Obviously, the 
perihelion shift can become quite large due to the conical nature of the space-time. Perihelion shifts in the solar system 
are quite small and the level of agreement of the measured value with that predicted by assuming the gravitational 
field of the sun to be described by the Schwarzschild solution can be used to give an upper bound on the energy 
density per unit length of a cosmic string present in the solar system. However, a large perihelion shift of 39 ° per 
orbit has been observed in a massive binary black hole system [19j| . Note that if we would do a similar approximation 
as that done in [ijj], we would find 

^=\{&4>)s + ^(\-l) , (38) 



where (A<f>)s denotes the perihelion shift in the Schwarzschild case (3 — 1. Within the PPN formalism [2(J, the 
perihelion shift can be written as Acj) = 42.98" Q(2 + 271 — 72)) assuming that the parameters 71 and 72 are equal 
to unity for General Relativity and that the quadrupole moment of the sun vanishes. Different experimental tests 
[HHU have given 72 - 1 = (1.2 ± 1.1) • 1(T 4 , hence s ^ 10~ 4 - 

If we now assume that the deviations of A0 from the Schwarzschild value are not due to a modification of General 
Relativity, but due to the presence of a cosmic string, we can approximate the energy per unit length of the cosmic 
string. Again reinstalling factors of c 2 we find (1 — /?) < 1CP 10 which for the deficit angle gives A < 10 -9 . This 
transfers to a bound on the energy per unit length which is given by: /j, < 10 17 — . 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

p 



FIG. 8: The perihelion shift in dependence on the deficit parameter /3 for m = 1, L z = 5, k = 1 and E — 0.999. 



C. Geodesies for D < and E 2 - |e > ^f-y 

This case corresponds to the choice V^ ax < ~ e (region III in the A-^t-plot). The effective potential and the 

corresponding characteristic polynomial P(z) are shown in Fig. [TJ] The line E ~ £ does not intersect the effective 
potential V s s(r). The polynomial P(z) has a single real zero at e\, however e\ < —1/3. We are allowed to integrate 
between z = 00 and z — —1/3, which corresponds to integration from r — to r = 00. 
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Effective potential Characteristic polynomial 




FIG. 9: The effective potential V e g{r) (left) and the characteristic polynomial P(z) (right) for D < and E 2 — |e > ■ 
The red horizontal line corresponds to the range of r values for an escape orbit. 



1. Escape orbit 

The range of allowed r values is indicated by the red horizontal line in Fig |5] and corresponds to < r < oo. The 
plot of an escape orbit in this case is shown in Fig llOl The particle comes in from infinity and - after a few loops - 
ends in r = 0. 




FIG. 10: We show a bound terminating orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin^, z — rcos8) 
(left), in the plane perpendicular to L (middle) and in R 3 (right). Here, we have chosen ft — 0.4, E = 2.8, L z — 5, k = 2.5 and 
m — 1. The blue circles in the p-z-plane and in the plane perpendicular to L correspond to circles with Schwarzschild radius 
r s — 2m = 2. 



D. Geodesies for D < and E 2 - |e < 

This case corresponds to the choice V^g m > — ~ e (region IV in the A-/i-plot). The effective potential and the 

corresponding characteristic polynomial P(z) are shown in Fig. 1111 The line E ~ £ intersects the effective potential 
V e g(r) only once at r\. This intersection point corresponds to the single real zero of the characteristic polynomial, 
which we denote by e\. Apparently, we are allowed to integrate between oo and e±, which corresponds to integration 
from r = to r = r\ . 
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Effective potential Characteristic polynomial 




FIG. 11: The effective potential V e s(r) (left) and the characteristic polynomial P(z) (right) for D < and E 2 — |e < ■ 
The red horizontal line corresponds to the range of r values for a bound terminating orbit. 



1. Bound terminating orbits 
The range of allowed r values is indicated by the red line in FigJTT] and is < r < r± . 




p 



FIG. 12: We show a bound terminating orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin^, z — rcosO) 
(left), in the plane perpendicular to L (middle) and in R 3 (right). Here, we have chosen /3 = 0.99, E = 0.9, L z = 5, 1/k = 0.99 
and m = 1. The dark yellow and blue circles in the p-z-plane and in the plane perpendicular to L correspond to circles with 
radius n and Schwarzschild radius r s = 2m = 2, respectively. 

The example of a bound terminating orbit for a massive test particle for j3 = 0.99 is given in Figfl^l 

E. Geodesies for D = and E 2 - e > 

This case corresponds to the choice V^ ax = ~ e > 0. The effective potential and the corresponding characteristic 

polynomial P(z) are shown in Fig. Q21 The line E 2 ~ E intersects the effective potential V e g(r) only once at n = r%, i.e. 
at the maximum of the effective potential. This intersection point corresponds to a double zero of the characteristic 
polynomial, which we denote by e\ — e-2.. Apparently, we are allowed to integrate between oo and e%, which corresponds 
to integration from r — to r = n and from e\ to z = —1/3, which corresponds to integration from r = r\ to r = oo. 
Note that the spherical orbit with r — r\ = ri = const, corresponds to the unstable spherical orbit discussed in 
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FIG. 13: The effective potential V e s(r) (left) and the characteristic polynomial P(z) (right) for D — and E 2 — e > 0. The 
red horizontal line corresponds to the range of r values for a bound terminating orbit, while the yellow line corresponds to the 
range of r for an unbound terminating orbit. 



1. Unbound terminating orbit 
The range of allowed r values is indicated by the yellow line in (|13[) and is r± < r < oo. 

An example of an unbound terminating orbit is given in Fig |14l for ft — 0.65. The test particle comes from infinity 
and approaches an unstable spherical orbit with r = r± = ?*2 = const.. 

2. Bound terminating orbit 

In this particular case the effective potential T4ff(r) has its maximum at the maximal radius r = r\ of the orbit. 
The range of allowed r values is indicated by the red line in (fT5|) and is < r < r\. 

An example of a bound terminating orbit is shown in Fig[l5] for j3 — 0.65. The particle starts at r — and 
approaches the unstable spherical orbit at r = r± = T2 = const.. 




FIG. 14: We show an unbound terminating orbit for a massive test particle (e = 1) in the p-z-plane (with p 
(left), in the plane perpendicular to L (middle) and in K 3 (right). We have chosen /3 = 0.65, E = 1.594, L z 
m = 1. The dark yellow circles in the p-z-plane and in the plane perpendicular to L correspond to circles 



= r sin 9, z = r cos 9) 
= 5, 1/jfc = 0.65 and 
with radius n = r<z. 
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FIG. 15: We show a bound terminating orbit for a massive test particle (e = 1) in the p-2-plane (with p = rsin6>, z — rcos8) 
(left), in the plane perpendicular to L (middle) and in R 3 (right). We have chosen j3 = 0.65, E = 1.594, L z = 5, 1/k = 0.65 
and m — 1. The dark yellow and blue circles in the p-z-plane and in the plane perpendicular to L correspond to circles with 
radius ri and Schwarzschild radius r s = 2m = 2, respectively. 



F. Geodesies for D = and E 2 - e < 

This case corresponds to the choice V^ ax = 2 ~ e < 0. The effective potential and the corresponding characteristic 

polynomial P{z) are shown in Fig. [TS] The line E ~ e intersects the effective potential V c g(r) at ri = r 2 and at 
r3. These intersection points correspond to a double zero ei = e 2 and a simple zero at 63. We are only allowed to 
integrate between e\ = e 2 and e3 which corresponds to integration between r\ — r 2 and 7^3 as well as between z = 00 
and ez which corresponds to integration between r = and T\ = r 2 . 

Effective potential Characteristic polynomial 




FIG. 16: The effective potential V e g(r) (left) and the characteristic polynomial P(z) (right) for D — and E 2 — e < 0. The 
red horizontal line corresponds to the range of r values for a bound terminating orbit, while the yellow line corresponds to the 
range of r for a bound orbit. 



1. Bound orbit 

The range of allowed r values is indicated in Fig ll6l bv the horizontal yellow line and is ri < r < r^. 
An example of a bound orbit is given in Fig ll7l for /3 = 0.9. This is an example of a so-called homoclinic orbit. 
Homoclinic orbits are orbits that approach unstable circular orbits and have been discussed extensively in the context 
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of gravitational wave production in binary system. There they play the role of a transition between the late inspiral 
and the final plunge |24| . 

2. Bound terminating orbit 

In this particular case the effective potential has its maximum at the maximal radius r = r\ of the orbit. An 
example of a bound terminating orbit is given in Fig ll8l 




FIG. 17: We show a bound orbit for a massive test particle (e = 1) in the p-z-plane (with p — rsinO, z = r cos#) (left), in the 
plane perpendicular to L (middle) and in K 3 (right). We have chosen /3 = 0.9, E — 0.9614, L z — 3.3, k = i2 anc J m = 1. The 
dark yellow and violet circles in the p-z-plane and in the plane perpendicular to L correspond to circles with radius r-2 and r%, 
respectively. 




FIG. 18: We show a bound terminating orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin^, z — rcos9) 
(left), in the plane perpendicular to L (middle) and in R 3 (right). We have chosen /3 = 0.9, E — 0.9614, L z = 3.3, k — ^ and 
m — 1. The dark yellow circles in the p-z-plane and in the plane perpendicular to L correspond to circles with radius ri. 
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G. Geodesies for D — and -Sr = ^} 

This case corresponds to the choice V^ m = — ~ £ . The effective potential and the corresponding characteristic 

polynomial P(z) are shown in Fig. Q1JJ The line E ~ £ intersects the effective potential KffM at r = 6m. This 
corresponds to e\ — e 2 — e 3 — 0. We are allowed to integrate from e\ = e 2 = e 3 to z = oo, which corresponds to 
integration from r± = r 2 = r% = 6m to r = 00. 

Effective potential Characteristic polynomial 




FIG. 19: The effective potential V e s(r) (left) and the characteristic polynomial P(z) (right) for D — and jj — t^^i- The 
red horizontal line corresponds to the range of r values for a bound terminating orbit. 



1. Bound terminating orbit 

The range of allowed r values is indicated by the horizontal red line in Figll9l and is < r < 6m. An example of a 
bound terminating orbit is given in Fig l20l for /3 = 0.9. 



2. Spherical orbit 

For ro = 6m we have a spherical orbit. These spherical orbits have already been discussed in [l6j . An example of 
a spherical orbit is given in FigJ21lfor j3 = 0.9. 



H. Geodesies for D = and E 2 - %e < -^-^ 

The effective potential and the corresponding characteristic polynomial P(z) are shown in Fig. [221 The line 
E ~ £ intersects the effective potential Kff( r ) at r\ and at the minimum of the effective potential at r = r 2 . These 
intersection points corresponds to a simple zero e\ and a double zero at e 2 = e^. We are only allowed to integrate 
between 00 and ei, which corresponds to integration from r = to r = r\. 

1. Bound terminating orbit 
An example of a bound terminating orbit is shown in Fig|23l 
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FIG. 20: We show a bound terminating orbit for a massive test particle (e = 1) in the p-2-plane (with p = rsin6>, z — rcos8) 
(left), in the plane perpendicular to L (middle) and in R 3 (right). We have chosen (3 = 0.9, E = a/8/9, L z = y/12 x (0.9) 2 , 
k = ^ anc ! m = 1. The dark yellow and blue circles in the p-z-plane and in the plane perpendicular to L correspond to circles 
with radius r = 6m = 6 and the Schwarzschild radius r s = 2m = 2, respectively. 




FIG. 21: We show a spherical orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin#, z — rcosO) (left), in 
the plane perpendicular to L (middle) and in R 3 (right). We have chosen j3 = 0.9, E — L z — a/12 x (0.9) 2 , k — ^ and 

m — 1. The dark yellow and blue circles in the p-2-plane and in the plane perpendicular to L correspond to circles with radius 
r = 6m = 6 and the Schwarzschild radius r s = 2m = 2, respectively. 



2. Spherical orbit 



Choosing ro = we find a spherical orbit with radius r2- In this case, the spherical orbit is a stable spherical orbit. 
The qualitative plot is similar to the plot given in FigETJ 



IV. NEWTONIAN LIMIT 



Assuming the gravitational field to be weak and the test particles to move slowly, we can give the Newtonian limit 
of our equations. This reads |26| : 

dr 2 1 - 

W = 2 Vh ™ (39) 
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FIG. 22: The effective potential V e g(r) (left) and the characteristic polynomial P(z) (right) for D = and E 2 — |e < ■ 
The red horizontal line corresponds to the range of r values for a bound terminating orbit. 





FIG. 23: We show a bound terminating orbit for a massive test particle (e = 1) in the p-z-plane (with p = rsin^, z — rcos9) 
(left), in the plane perpendicular to L (middle) and in M 3 (right). We have chosen j3 = 0.9, E — 0.95229, L z = 3.3 and k = 



where we assume g^u — rj^ u + with very small and i]^ to be the Minkowski metric. Proceeding as in classical 
mechanics, we find the Virial theorem 



2 < T >= i < r- V^oo >=\< rd rhoo > 

2 



(40) 



where < .. > denotes the temporal average and T = | {%:)* is the kinetic energy. Hence, the Virial theorem doesn't 
change for (3 ^ 1 as compared to the Newtonian limit of the standard Schwarzschild case (3 = 1. However, note that 
the kinetic energy T has a /3-dependence since 



r £ sin 2 Of 



dt 



(41) 



Since the average of the potential energy < rd r hoo > does not depend on /3 and the Virial theorem holds for all /3 this 
means that when /3 < 1 the angular velocity {^^j increases. We have - in fact - seen this effect when studying the 

perihelion shift of planets. For the same interval of eigentime, the planet moves further in a space-time with j3 < 1 as 
compared to in a space-time with = 1. 

In addition, it is easy to show that Kepler's third law T 2 /a 3 = const., where T is the period of the orbit and a the 
semimajor axis also doesn't change. This is easiest to see when considering circular and planar orbits with L z = L, 
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i.e. — 7r/2. For these, we let r = and substitute L z /r from (j> into r (see (fTTj) and ([12]) ). After separation the two 
sides can be integrated to give the period of the orbit. The relation will contain a dependence on /3. However, note 
that the integration in 4> is from to 2tt(3 such that the /3 cancels and T and a are independent of /?. 

V. CONCLUSIONS AND OUTLOOK 

In this paper, we have constructed the complete set of geodesies in the space-time of a Schwarzschild black hole 
pierced by a cosmic string. We have classified the geodesies according to their energy, angular momentum aligned 
with the axis of the string L z , the ratio between total angular momentum and L z , the horizon radius of the black 
hole and the deficit angle (or energy per unit length of the string). We found that the change of the deficit angle 
doesn't effect the r(ff) motion, but that the r(0) motion changes significantly. Moreover, the motion of the particles 
is in general not planar. We found a bound of the energy per unit length of the string comparing our results 
with experimental tests of General Relativity. From perihelion shift and light deflection we found fi < 10 17 ^ and 

M < 10 16 ^S, respectively. The existence of cosmic strings may also influence the creation of gravitational waves since 
for inspirals and flybys the deficit angle modifies the angular velocity and, thus, the temporal change of the effective 
quadrupole responsible for the emitted gravitational wave. The results obtained in this paper can be generalized 
to rotating black holes, i.e. to the space-time of a Kerr black hole pierced by a cosmic string. They can further be 
generalized to a Schwarzschild-de Sitter black hole pierced by a cosmic string. For treating such orbits one has to 
apply hyperelliptic integration developed in ||. We expect the same enlarged variety of orbits as we encountered in 
Schwarzschild-de Sitter space-time as compared to Schwarzschild space-time. 
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